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■ We derive the initial condition for the perturbative fragmentation function of a heavy quark 
' through order 0{a 2 s ) in the MS scheme. This initial condition is useful for computing heavy quark 
\ (or lepton, in case of QED) energy distributions from calculations with massless partons. In addition, 

■ the initial condition at 0{a 2 s ) can be used to resum collinear logarithms ln(Q 2 /m 2 ) in heavy quark 
(-H ■ energy spectrum with next-to-next-to-leading logarithmic accuracy by solving the DGLAP equation. 

>—i . 

(T) ; I. INTRODUCTION 

i ! 

1 | Production of heavy flavors (charm and bottom) in high energy processes has become an important subject in the 

■ last decade, due to experiments at e + e~ and hadron colliders. Since the center of mass energy of these machines is much 
£f*l | larger than bottom and charm quark masses, it is tempting to consider such processes in the massless approximation. 
t^J- , Unfortunately, this is only possible for sufficiently inclusive observables, for example, total cross-sections. However, 
t-H ' differential distributions are more relevant for experimental analysis and, in addition, contain more information about 

the underlying physics. An interesting example is the energy distribution of heavy hadrons produced in a high-energy 
collision, since it gives direct access to the hadronization process. In perturbation theory, however, the heavy quark 
energy distribution diverges in the limit to — > and the sensitivity to the quark mass m remains at arbitrary high 
energies. 

This sensitivity is disagreeable for two reasons. First, it implies that energy spectra must be computed retaining 
quark masses which makes the calculations quite involved. Second, in higher orders of perturbation theory, energy 
O ' spectra contain powers of logarithms ln(Q 2 /m 2 ), that may invalidate fixed order perturbative calculations. 

, Both of these problems are solved by introducing perturbative fragmentation function [| . This function describes 
' the probability that a massless parton of a certain type fragments into a heavy quark. The mass of the heavy quark 
provides a natural cut-off for the collinear radiation, thus making perturbative predictions finite. The perturbative 
fragmentation function satisfies the DGLAP evolution equation |2(. By solving this equation, we can sum up large 
collinear logarithms ln(Q 2 /m 2 ) that appear in perturbative fragmentation function and improve the quality of pertur- 
bative calculations. To avoid confusion, we note that a meaningful prediction for heavy hadron energy spectrum can 
only be obtained if the heavy quark energy distribution is convoluted with a non-pcrturbativc fragmentation function. 
Such functions are traditionally extracted by fitting hadron energy spectra in e + e _ collisions. We will not discuss 
this issue here; for recent work on the subject see j3|- We also note that the concept of perturbative fragmentation is 
useful outside QCD; for example, QED effects in the electron energy spectrum in muon decay, currently under study 
by the TWIST collaboration £| , can be analyzed using the same technique |f| . 

Because the Q 2 evolution of the fragmentation function can be obtained from the DGLAP equation, the perturbative 
fragmentation function at arbitrary Q 2 can be fully reconstructed once the initial condition D ml and the time-like 
splitting kernels for the DGLAP evolution are known. The initial condition is usually computed at Q 2 ~ to 2 ; this 
ensures that no large logarithms appear and fixed order perturbation theory is reliable. 

The purpose of this paper is to provide such initial condition for the perturbative fragmentation function through 
0(a 2 ). As we explained earlier, the knowledge of D lm through next-to- next-to-leading order (NNLO) enables the 
calculation of energy spectra of heavy quarks through 0{a 2 ) from massless results and, simultaneously, allows to 
resum large logarithms ln(Q 2 /m 2 ) using the DGLAP equation. All power corrections 0{m 2 /Q 2 ) are neglected within 
this approach, but for many practical purposes such precision is sufficient. 

Technically, the use of perturbative fragmentation functions makes it possible to compute the heavy quark energy 
spectrum in two steps - we first compute the energy spectrum of massless partons produced in a hard scattering and 
then convolute this spectrum with perturbative fragmentation function. The possibility to remove all the dependence 
on the heavy quark mass from the calculation of the hard scattering simplifies the computations considerably. 
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The remainder of the paper is organized as follows. In Section II we describe the process-independent derivation 
of the initial condition for the fragmentation function. Then, we discuss collinear factorization when finite quark 
masses are present and explain how to compute D lm from relevant Feynman diagrams. We also suggest a suitable 
modification of the original proposal [(J 0] for the process-independent computation of the initial condition; such 
modification significantly simplifies NNLO calculations. In Section IV the initial condition for the fragmentation 
function for both quark- and gluon-initiated processes is computed through NLO. This computation allows us to 
demonstrate the details of our approach and to derive the NLO perturbative fragmentation function through O(e), 
where e is the dimensional regularization parameter. In Section V we describe the calculation of the 0(a 2 s ) contribution 
to the fragmentation function and present the result for the initial condition. We conclude in Section VI. 



II. PROCESS-INDEPENDENT DERIVATION OF D ini 



Consider production of a heavy quark Q with mass m and a definite value of energy Eq in a hard scattering 
process. According to the QCD factorization theorems |sll9UlC|. the heavy quark energy spectrum can be computed 
as a convolution of the energy distribution of massless partons produced in the hard process, and the fragmentation 
function that describes the probability that the massless parton fragments into a massive quark with a definite energy. 
If the energy fraction Eq/ EQ_ max of the heavy quark is denoted by z, then the energy distribution of that quark can 
be written as: 



dcr Q 

dz 



E 



1 dx do, 

x dx 



-(x,Q,n)D a/Q (-, — ) 

C \X 771/ 



(1) 



Here the sum runs over all partons (quarks, antiquarks and gluons) that can be produced in the hard process and p, is 
the factorization scale. The coefficient function da a jdx is the MS renormalized differential cross-section for producing 
a massless parton a 1 . It is defined indirectly through the equation 



d<T a 

dz 



(z, Q,e) = Y^ ~p-( z i Q> M) ® r ba (z, fi, e), 



(2) 



where da a /dz is the bare energy distribution for the parton of type a; the collinear divergences in this distribution 
are regularized by working in d = 4 — 2e dimensions. T a i, are universal collinear subtraction terms, defined in the MS 
scheme: 



Tfca = S ab S(l - z) 



fOLs_ 
V27T 



V2vr 



(3) 



where a s = a s (/z) is the MS strong coupling constant, renormalized at the scale /i. The relation between the bare 
and the renormalized couplings reads: 



2tt 



a s 
2^ 



2ttT 



0{a 2 s 



(4) 



where S £ = (Air)' 
/3-function, Ca z 



e 67 and 7 is the Euler constant. Also, (3q = (11CU — 4T^n/)/6 is the 0(a 2 s ) coefficient of the QCD 
= 3, Tr = 1/2 are the QCD color factors, n/ denotes the number of fermion flavors (including Q) 



and -F^b ^ are ^ ne time- like splitting functions [Tl|. Our notations for the splitting functions follow Ref.[T(i| 

The functions D a /Q(x, fi/m) in Eq.lJIJ are the perturbative fragmentation functions 0. They satisfy the DGLAP 



evolution equation and can be fully reconstructed from it, if the initial condition at a scale fi = /io is known, 
denote 



D 



a/Q 



V ' m / a V ' m / 



We 



(5) 



In the evaluation of the coefficient function da a /dx the heavy quark Q is considered as massless; therefore, the sum over indexes in 
Eq. includes the flavor Q. 
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If Hq ~ m is chosen, the initial condition D™ 1 does not contain large logarithms and can be derived from fixed order 
perturbative calculations. Setting fi = fi in Eq.JQ), we obtain: 



da 



dz 



Q,m)=Y J ^ri*, Q, Mo) ® £>f (*, ^) 
j * — ' dz V m / 



(6) 



where dcr fo is the fixed order cross-section for producing a quark Q with mass m, with all power corrections 
0((m 2 /Q 2 ) n ), n > 1 systematically neglected. It follows from this equation that the initial condition for the frag- 
mentation function can be extracted, if the heavy quark energy spectrum is computed in a particular process through 
certain order in a s . This, however, is technically inconvenient. 

A more convenient, process-independent approach was suggested recently in Refs.0, Q- Since hard scattering 
cross-sections are insensitive to long-distance dynamics, we may write 



—^(z,Q,m) = 2^ -T- {z,Q,no) <g> D a/Q 



dz 



Mo 



— {z, Q,e) = -1-{ Z > Q> Mo) ® D a/b {z, Ho,e). 



(7) 
(8) 



The functions D and D L describe collinear radiation from massive (massless) quarks; we will show below that they 
can be considered as bare fragmentation functions. Combining Eas, 17l8|) with Eas. i|6l2|) . we obtain: 



Da/Q (z, ^) = E 5 a/b(z, MO, e) ® [T( 
b,c 



z.Mo,e)J 6c 



■m 



(zA) 



(9) 



It follows from this equation that D ml can be derived once D and D L are known. We describe how to compute D 
and D L in the next Section. 

To conclude this Section, we comment on the analytic structure of D lnl . The initial condition for the perturbative 
fragmentation function can be expanded in series of a s : 



n=0 v 7 



The results through 0(a s ) are known [lj: 

d£ 0) (*) = 



1 + z 2 



1 - z 



111 



Mo 



4 1 2 s (.^)=T fl (z 2 + (l-z) 2 )ln(^ 



m 2 (l — z) 1 



Mo 



fz ^ 
d a7tQj9 [z, I 



= 0, 



(10) 



(11) 



where Cf ~ 4/3 is the QCD color factor. In this paper, we are mostly concerned with the computation of the 
coefficients df® — _, where q stands for any massless quark. The logarithmically enhanced terms in d^ follow from 
the DGLAP equation: 



d{ 



i 2) (*, 



Mo\ 
m ) 



P 



(0) 



ba 



l P Qb( Z ) A)p(O)/ x 

~2 + ~2 Qa( ' 



hi" 



Mo 



Pali*) + ® 4 X) (*. 1) + ft^fr 1)1 In ( Q ) + d< 2) (*, 1) 



(12) 



The mass-independent constant of integration da (z, 1) can not be determined from the DGLAP equation and has to 
be computed explicitly. This is the major goal of this paper. In the next Section we describe how this can be done. 
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III. THE COLLINEAR LIMIT 

The idea that allows explicit computation of the functions D L and D and, therefore, of the initial condition for 
the perturbative fragmentation function is as follows. Both D and D L are introduced to describe collinear radiation 
with q± ~ to < fj,Q- Since we are not interested in power-suppressed contributions to the cross-section, we can restrict 
ourselves to such terms in scattering amplitudes that, upon integration over the phase-space of final state particles, 
produce C(ln"(TO)), n > 1 terms 2 . Such terms can be identified by applying power counting arguments [3 in the 
collinear limit for scattering amplitudes. The collinear limit is defined as the kinematic limit where the relative 
transverse momentum of two or more particles vanishes. When masses are introduced, the collinear limit is q± — > 0, 
to — > and q±/m = const [T^ . 

Throughout the paper we work in the light-cone gauge n^A^ = 0, n 2 = 0. This is necessary for the process- 
independent derivation of D lni because in such gauges, the 0(ln n (?7i)), n > 1 terms are produced only by diagrams 
where collinear radiation is both emitted and absorbed by the same parton. The quantum interference effects in such 
gauges can be integrated over the phase-space of the final state particles in the massless approximation. Because of 
that, D lm can be computed from self-energy type diagrams, integrated over the virtuality of the incoming parton. 
Different cuts of such self-energy diagrams correspond to different contributions to the fragmentation function: at 
0(a 2 s ), we have to deal with two-loop virtual corrections, one-loop virtual corrections to one-to-two splittings and, 
finally, one-to-three splittings. 

We begin by considering the one-to-three collinear splitting since the kinematics of the final state in this case is the 
most general. We denote the four-momenta of the produced (massive or massless) particles by 91,2,3 and their sum 
by p : 

p = qi+q 2 + q 3 - (13) 

We parameterize the collinear direction by p and use the gauge fixing vector n as the complimentary light-cone vector. 
We then write: 

Qi = ZiP + fc- — r+qi,±. (14) 
(pn) 

The components f3i are found from the on-shcll conditions q 2 — m 2 : 

-q? 1 + mf - z 2 to 2 

A= — 9 — — , »= 1,2,3, (15) 

with p 2 = to 2 . We are interested in the case when p, p and n belong to the same hyperplane in the ci-dimensional 
space 3 . Then: 

zx + Z2 + z 3 = 1 , qi t ± + q 2 ,x + Q3,± = 0. (16) 
Using the above equations, we express the momentum p through p and n: 

where: 

^2 2 2 (?"?2) + 2(pq 3 ) - 2{q 2 q 3 ) + m\ - m\ - m\ - m 2 
P ~m = ^. 

' z 

Our notations are such that qi always denotes the momentum of the heavy quark in the final state, whose energy is 
measured. For that reason we often write z instead of Z\, to denote its energy fraction. Eas. (|f 411 6|) imply the relation: 

, , (nga) (nq 3 ) 

z - 1 - 7 — \ 7 — r, U 8 ; 

(pn) (pn) 



2 To be precise, we are interested in all contributions to the scattering amplitude that, upon integration over the phase space, scale as 
m~ nE , where n is some integer. Upon expanding in e, those terms produce both ln(m)-enhanced and m-independent contributions to 
the energy spectrum. 

3 More general configuration has been considered in ll.'l . 
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which is nothing but the energy conservation condition. 

Although we have only discussed the kinematics of the one-to-three splitting, the kinematics of the one-to-two 
splitting can be obtained as a particular case. For this, it is sufficient to set the momentum q 3 and the mass TO3 to 
zero in the equations above. 

In QCD, factorization properties of the cross-sections, as in Eas. (|7l8[) . can be traced back to the factorization 
properties of both the phase-space and the matrix elements in the collinear kinematics. We illustrate those properties 
by considering the tree level amplitudes and then generalize these considerations to include virtual corrections. 

Consider a hard scattering process characterized by some scale Q 3> to. Suppose that n + 2 partons with momenta 
k%, . . . , fc„_i, qi, q 2 , 93 are produced. We consider momenta k\,..., fe n -i as non-exceptional, while momenta qi, (72, 53 
are collinear, as described above. First, we discuss the phase-space factorization. 

We denote the phase-space of n + 2 partons as dPS^™ +2 ^(fci, . . . , fc n _i, q%, q 2 , 93)- In the limit when qi, q 2 and 93 
become collinear qi + qi + q 3 = p + 0(q±), the (n + 2)-particle phase space factorizes: 

dPS (n+2) (fc!, . . . , fc„_i, 9l, 92, 93) = dPS< n >(fci, • • • , K-X,p) d<S> col \q 2 , q 3 ). (19) 
We use the following notations: 

d$ coll (<Z2,g 3 ) = -[dq2;m 2 ][dq 3 ;m 3 }, (20) 

z 

and [dq; m q ] is the (i-dimensional one-particle phase space: 

[dq;m q }^j0zrS + (q 2 -m 2 q ). (21) 

To derive Eq.lplJ), we neglect the 0(q±,m) difference between p and p in dPS (n) and integrate over qi using Ea. (|13f) . 

We now consider factorization properties of the tree-level matrix element. In the massless case, they were described 
in |l3lll4l]. Below we generalize those results to the case of non-zero mass m. We consider massive quarks produced 
in the fragmentation of a quark-like parton of either the same or different flavor. The case of a gluon fragmentation 
into a heavy quark can be dealt with in a similar way. 

We already mentioned that, in physical gauges, the relevant contributions in the collinear limit come from di- 
agrams where the same parton emits and absorbs collinear radiation. Taken in conjunction with power counting 
arguments, this observation can be used to simplify tree-level matrix elements in the collinear limit and derive the 
factorization formula. Consider a process with n + 2 particles in the final state described by the matrix element 
M^ l+2 \ki, . . . , k n -i, qi, (72, 93)- In physical gauges, the collinear splitting p — > q\ + q 2 +03 decouples from the rest of 
the process, so that the squared, spin-averaged amplitude |Af(™ +2 )| 2 can be written as |13| : 

\M^ +2 \k x , kn- u 91, g 2 , q 3 )\ 2 =M^(k u ..., k n - 1 ,p) fl , a VZ$(p, 92, 93 ). (22) 

Here a, (3 are spinor indices. If, in the collinear limit, V c ° n can be written as 

V*f =i a3 W + 0{ qi _), (23) 

we can rewrite Ea. (|22|l to make the factorization of the matrix clement explicit: 

iM^+^fcx, . . . , fe n _!, q x , q 2 , q 3 )\ 2 = \M (n \k u fc„_i,p)| 2 W(n,p, q 2 ,q 3 ) + 0(q±). (24) 

Here | (...., p)\ 2 is the amplitude squared for producing n on-shell particles with non-exceptional momenta; the 
scalar function W contains all the information about the collinear splitting. 

We now explain why V c ° n can be written as in Eq. (|23|l . Because the collinear splitting p — > q\ + q 2 + q 3 decouples 
from the rest of the process, V c ° n is obtained as the sum squared of all possible diagrams that describe a transition of 
an off-shell parton with momentum p into three particles with momenta 91,2,3- The propagator of the off-shell parton 
p is included in V^ co11 . In the massless case the matrix V coil (n,p, q 2 , 93) can be written as jl3j : 

v-"SfHf?'* +B (^ 2 )- <25) 

The functions Aj and B are dimensionless functions that depend on the scalar products (9^9^). To determine the 
behavior of the amplitude V co11 in the collinear limit, we rescale the transverse momenta: 



q it ± -> Kq it ± (26) 
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and consider the limit k — > 0. Since (qiqj) ~ O(g^), the functions Ai_ 3 and -B remain invariant under that rescaling. 
As follows from Ea. l|14l) . the matrices 4i transform as: d { — > z, fi + O(k). Therefore, after the rescaling Eci. (|26[l . the 
amplitude V c ° n becomes: 

r'"" = „ '['-^) I >M,:, I P C)(/, i .27) 




Because the collinear phase space Ea. (|2U|) transforms as: 

d$ coll (g 2 , g 3 ) - n 4 -^d^ oU (q 2 ,q 3 ), (28) 

under the rescaling Eq. l|26l) . it follows that in the collinear limit k — > only the first term in Ea. (|27|l gives non- 
vanishing contribution after the phase-space integration. This proves that V co11 has the form shown in Ea. (|23[l . We 
can extract the function W from V c ° n by applying the projection: 

Tr \riV coU ] 

We now generalize this result to the case when at least one of the final state particles is a quark with mass to. The 
mass to is considered to be of the order of the transverse momenta of the collinear particles. We follow the same line 
of reasoning that leads to Ea. (|25ll allowing, however, for non- vanishing masses in the initial and final states of the 
collinear splitting process. Ea. (|25ll generalizes to: 

In Ea. l|3U|) . rrii n is the mass of the parton that initiates the collinear splitting. The functions A[ m \ £?( m ) and C 
are scalar functions of dimension zero, two and zero respectively; they depend on (qiqj) and to 2 . These functions 
possess a regular to — > limit. In the collinear limit, when both the transverse momenta and the mass to are rescaled 
simultaneously, 

9i,_L ► K 9i,X) m — > k to, (31) 

the scaling properties of the functions A[ m \ _B( m ) and C coincide with their mass dimension. Therefore, under the 
rescaling Ea. (|31|l . the amplitude y co11 behaves as: 

2£ \ 2 / 3 \ 

3 )- (32) 

This result implies that also for the to ^ case, the function W is given by Ea. (|29|) . Finally, assembling all the pieces, 
we find that the function D(z) in Eq.JJJ), can be written as: 

D(Z) = - J [dq 2 ; m2 ][dq 3 ; m3 ] ^ 5 \l z — — j . (33) 

Although we have derived Ea. (|33(l for the one-to-three splitting contribution to D™ 1 (z, fXo/m), similar expression is 
valid when virtual corrections are included. In addition, Eq. I|33J) can be used to derive the function D L after the mass 
to is set to zero everywhere. 

We now discuss how to perform phase space integrations in Eq. l|33|l . The approach proposed in 0, Q requires 
integration over the transverse momentum up to the scale = /io for both D and D L . While this allows a simple 
calculation of the fragmentation function at 0(a s ), this approach becomes impractical in higher orders of perturbation 
theory. It is possible to simplify the calculation considerably by realizing that the limit fi± — > oo can be taken when 
the difference of massive and massless functions D is considered. The technical simplification associated with this 
is twofold. First, all higher order QCD effects in the massless function D L vanish, because all the integrals that 
contribute there become scaleless. We obtain: 

D^ /b (z,fx ,e)=S ab S(l-z). (34) 




I I 

Q ' Q I 



FIG. 1: Diagrams that contribute to perturbative fragmentation of a heavy quark Q — > Q + X at 0(a s ). The dashed vertical 
line indicates the intermediate state that has to be considered. 

Second, in the massive case, only single-scale integrals have to be computed. Using Ea. l|34fl . we obtain a simple 
expression for the initial condition of the perturbative fragmentation function: 



DT U —) = E r ^> «>) ® Db/Q U ^) . (35) 
V ml *■ — ' V ml 



m I 

b 

From this it follows, that D b /Q is the bare fragmentation function for the massive quark, whereas D lnl represents its 
collinear renormalized version. 

As we mentioned earlier, there are three contributions to D b /Q(z, to), that have to be considered: double- virtual, 
real-virtual and double-real. In all these cases the amplitude V c ° n can be constructed by applying conventional 
Feynman rules to describe the transition of a parton with momentum p to the final state of up to three particles 
with momenta 91,2,3- Virtual corrections should also be included, when appropriate. For each of these three cases, 
the function W is obtained using the same projector as in Ea. l|35f) . The only unusual feature of the amplitudes that 
contribute to y co11 is that the spinor, that describes the initial state parton, has to be replaced with the propagator 
of the same parton with the off-shell momentum p. V co11 is then given by the square of the corresponding matrix 
element. Standard symmetry factors apply in this case as well. In particular, to evaluate D one has to sum over 
spins and colors of the final state and to average over colors in the initial state. The ultraviolet renormalization is 
performed in a standard fashion. One has to renormalize the QCD coupling constant, the quark mass and external 
quark and gluon fields. Since we work in the light cone gauge, we find it convenient to compute the quark and gluon 
field renormalization constants by computing diagrams with self-energy insertions to external quark and gluon lines. 

Finally, we comment on the double virtual contributions to the fragmentation function. Such contributions describe 
the one-to-one splitting and are proportional to 5(1— z). Unfortunately, computing these contributions requires dealing 
with well-known complications of light-cone gauges. For this reason, we decided to fix the (e-dependent) coefficient 
of 5(1 — z) using the fermion number conservation sum rule: 

dz{D%) Q (z)-Df /Q )=l. (36) 







We do not encounter, however, any spurious light-cone singularities in the diagrams with one virtual loop. All 
singularities originating from those diagrams are consistently regulated with dimensional regularization. 

To perform integrations in phase-space and loop integrals we proceed as in |l5l Il6l Il7| , where further details of 
the method can be found. The idea is to map all non-trivial phase-space integrals to loop integrals and use standard 
multiloop methods [Tsj . such as integration-by-parts and recurrence relations, to reduce all the phase-space and loop 
integrals that have to be evaluated, to a few master integrals. For the reduction to master integrals we use the 
algorithm j^j implemented in [2(j . All algebraic manipulations have been performed using Maple [2JJ and Form • 
We now discuss the calculation of the initial condition D lnl through 0(a s ) which allows us to illustrate the details of 
the method by considering a simple example. 

IV. D ini AT NLO 

In this Section we compute the initial condition for the fragmentation function at NLO through 0(e). Such terms 
are needed for the computation of D lm at order 0(a 2 s ). The function V coU for Q — > Q + g splitting is obtained by 
considering the diagrams shown in Fig.l. To facilitate the comparison with the literature, we express the function W, 
Ea.p§)l. through the so-called splitting function P [l3[ : 

\£r — m ) 

Introducing auxiliary vector n: 
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and using the projector Ea. (|29() . we obtain the splitting function: 



(qn) 



2 — m 



2 (pn) - (w) 



(pn) ' (pq)(pn) 

Combining Eas. (|18l38|l . we derive the real emission contribution to the fragmentation function D: 



2, J 

Because of the constraint 5 (1 — (nq)) in Ea. (|40|l . the splitting function simplifies: 

2 



1 - z 



— 1 — z — e(l — z) — z 



(pq) 



(39) 



(40) 



(41) 



The result for pW in Eq.^ agrees with Ref.[T3|- When P^ is used in Ea. (j4*U|) we observe that two integrals have 
to be evaluated. However, those integrals are not independent; an algebraic relation between them can be found using 
the method of 

HEHEQ As a result, D^(z) can be expressed through a single "master" integral: 



We finally obtain the real emission contribution to the fragmentation function: 



V mJ \2ir J \ m 2 J (1 - z) 1+2e 



valid to all orders in e. To expand this result in e, we use 



(1-z)- 



E- 

n>0 



ln"(l -z) 
1 - z 



(42) 



(43) 



(44) 



Virtual corrections are derived by considering self-energy diagrams in the light-cone gauge. The result reads: 



Mo\ 
m ) 



v (4 



2tt 
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1 1 



S e C F 6(l-z)<- 5 + 



Mo 



In 



Mo 



Vf4 

m 



Mo 



2 

C(3) 



12 



+ £>(e 2 ). 



6 \m 2 J 4 Vm 2 y V 12 / \m 2 J 24 3 
Combining Eas. (|43l45|l . we arrive at the bare perturbative fragmentation function through 0(a s 



(45) 



1 - z 



■In" 



Mo 



-In 



Mo 



1 + z 2 



Mo 



In inL -21n(l-z)-l 



(2 ln(l - z) + 1) + 2 ln 2 (l - z) + 2 ln(l - z) + 



12 



0(e 2 



(46) 



As expected from the fermion number conservation condition, Ea. (|36|l . the integral of Dq/q(z, fio/m) over z equals 
1 and does not receive any corrections at order 0(a s ). Applying collinear renormalization Ea. l35|) and substituting 
a®S e — > a s in Ea. (|46[) . we observe that the 1/e term in Ea. (|46|) cancels and the remaining terms are not modified. 
Therefore D lm can be read off from Ea. (|46[) by simply omitting the 1/e term. The O(e ) part of D ml coincides with 
the result of Ref. |3, while the 0(e) part is new. 

The gluon fragmentation g — > Q + Q can be treated in a similar way. The real emission contribution reads: 



Q Q < Q 

■ i i 



FIG. 2: Examples of diagrams that contribute to perturbative fragmentation of a heavy quark Q — » Q + X at C?(a;?). The 
dashed vertical line indicates the intermediate state that has to be considered. 



Evaluation of the splitting function P g /Q is described in [f|. We derive: 

„ „ / 1 — z z m 

which agrees with the result in Ref. There are two integrals of the type: 

J(a) = / w J w - m2] 5[1 - i?iq)] ' a=1 ' 2 

that have to be computed, but only one of them is independent. We obtain: 

J (1) = ( 7 r) 1 -T(e)m- 2£ (l - z). 



(48) 



(49) 



(50) 



This allows us to derive the gluon-initiated contribution to the initial condition for the perturbative fragmentation 
function at 0(a s ): 



lln 2 (^) + - 



0(e 2 



(51) 



The 1/e pole cancels after collinear renormalization and the known result for D g /Q Q is recovered in the limit e — > 0. 
Having explained out method by considering a simple example, we now present the results for D lnl through 0{a 2 s ). 



V. D ini AT NNLO 



At NNLO, large number of processes contribute, yet, all of them can be treated using the methods described above. 
Before we present the results for D lnl through 0(a 2 ), we briefly mention some technical details that we find peculiar. 

Several sub-processes contribute at 0(a 2 s ). It is convenient to split them according to their contributions to various 
components of D 1 ™. The real emission contributions to Dq 1 are Q — ► Qgg. Qqq, QQQ, while for D^ 1 and D™^ the 

real emission sub-processes are Q — > QQQ and q(q) — > Qq(q)Q, respectively. Some of the contributing diagrams are 
shown on Fig. 2. 

Because the number of contributing processes is large, there are many possibilities of mass assignments for particles 
in the final state, relevant for the calculation of D lnl . We have to consider all of those cases separately. As we 
explained in Section III, we require the tree-level splitting amplitudes in the collinear limit when massive quarks are 
present. Those results are not available in the literature; yet, a useful cross-check of our results is obtained once the 
limit m — > is taken. In that limit the splitting amplitudes derived in this paper coincide with those in Ref.|l3|. 

We now present the final results for the fermion initiated contributions to the initial condition of the perturbative 
fragmentation function at order a 2 .. We give the results for the coefficients (z, fio/m), introduced in Ea. (|10|l . Our 
results contain polylogarithmic functions up to rank three. These functions are defined through: 

Lin (2) = / — — ] —-dx, Lii(z) = — ln(l — z), Si, 2(2) = - / — -dx. (52) 

Jo x 2 J x 

We begin with the simplest case of a light (anti)quark-initiated fragmentation process. Through 0(a 2 ) : the two 
contributions coincide: 

(••£)■ < 53 > 
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We write: 

d^=C F T R Ff- T -\ (54) 

where 

F?™ = {(l + S )ln W+ "-" 4 ; 2 3 + " + 4 » } L> + {(1 + z) ln*<s) - fc2 + 2 ^ + ' 5 ln( 2 ) 
_ 4 (1 -„ (1 ^ + 2 3 , + 5 ) | L + ^ + 2) + (l-,)(^ + 7, + 4) ) (L . (>) + ¥1 _ z) 



+ 4(1 + 2) 



Li 3 (z) - ln(z)Li 2 (z) + 1 ln 3 (z) - 1 ln(l - z) ln 2 (z) - C(3) 



16 , 56_ 56^,^^ , 7(1 - z)(16 + 133z + 88z 2 ) 



15 + 27z + 8z 2 2 
12 ln W 

+ ( f + ^ + ^ J Hz) + - ^- — ^ (55) 
where L — ln^/m 2 ). 



Next, we present the heavy antiquark-initiated contribution c0 . We write: 



4j } = c F (c F £±) + ^4^. 



The functions f!£ a,Cf) and Fi CFTi,) read: 

4 Ca ' Cf) = (-4Li 2 (-z) + ln 2 (z)-41n(l + z)ln(z)-y) + 2(1 + z) ln(z) + 4(1 - z) ^ L 

1 + z 2 / \n(z) 3 
+ - 1 — 4Li 3 (z) - 2Li 3 (-z) + 12S h2 (-z) + ~Y~ - 2 ln(z)Li 2 (z) - 2 ln(z)Li 2 (-z) 

27T 2 

+ 12 ln(l + z)Li 2 (-z) - 3 ln(l + z) ln(z) 2 + 6 ln 2 (l + z) ln(z) - — ln(z) - 7C(3) + tt 2 ln(l + z) 
4(12z 2 + 3z 4 + 8z 3 + l + 12z) T . . . 2(3z 2 - 1 - 10z) T . . . z(-3z + 3 + z 2 + 3z 3 ) , ,. , 

3(TT^ Ll2(z) + 3(1 + z) Ll2( - z) 3(TT^ ln {z) 

4(12z 2 + 3z 4 + 8z 3 + l + 12z) 1 2(3z 2 - 1 - lOz) 

^TT^ ln d - *) + — 3(iT7) — ln(1 + z) Hz) 

-30z 2 + 15z 4 - 9 - 28z + 36z 3 30z 2 + 28z 3 + 15z 4 + 36z + 3 2 45z 3 - 29z + 29z 2 - 45 

+ 6(l + z)3 {Z)+ 18(1 + z) 3 * + 6(1 + z) 2 ' 



and 



(56) 



(57) 



r,(c F T R ) („ NW , 4z 2 + 3z-3 2 1 2 f 2 8z 2 + 27z+15, , . 
4 = [(l + ^)ln(z) + -|L 2 + |(l + z)ln 2 (z) ln(z) 

+ 56z 2 + 36z-72 _ 20 j £ + (1 + ^ ^ 41n(z)Li2(z) + 81n(z)Li2 (_ z) _ 8Li3(z) _ i 6 Li 3 (-z) + iln 3 (z) - 4C(3)) 

+ 4(Z ~ 1)(Z 2 + 4Z + 1) (Li 2 (z) + 2Li 2 (-z) + ln(l - z) ln(z) + 2 ln(l + z) ln(z)) - gf!±j!i±^ ln 2 (z) 
oz 12 

20z 5 + 351z + 48 + 147z 4 + 489z 3 + 713z 2 , , , 400z 4 - 963z + 1331z 3 - 819 + 499z 2 56 

+ 9(IT7)3 ln(z) mTzf + 27? (58) 



11 



Finally, we present the heavy quark-initiated contribution cIq\z, /Uo/m). We write: 



4 2) U- 



C 2 F F^ l) + C a C f F^ aCf) + C f T r F<£ fTr) + C F T Rni F^ FTRni \ 



where ri; = nj — 1 is the number of massless flavors. The functions Fq read: 



(59) 



p(C|) _ 



8 3 



(5(1 -z) + 4 



ln(l - z) 



1 - z 



5 + z 



27 7T 



L + {[ "5- + ir - 2 C(3) <5(1 - ^) - 12 



+ 3 
- + 


1 


1 - z_ 


- 12 


'ln 2 (l - z 


1 - z 



l + 3z 2 



J + 



2(1 -z) 

"ln(l - z) 



ln(z)-2(l + z)ln(l-z) 
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1 - z 



+ 1 + 



4tt' 



1 - z 



+ (1 + z) (2Li 2 (z) -7T 2 ) 



3 + 5z 2 
2(1 - z) 



ln 2 (z) 



+6(1 + z) ln 2 (l - z) + (11 + 3z) ln(l - z) + 



6(1 + z 2 ) 
(1-z) 
9z - 11 



ln(l - z)ln(z) 



(4z 2 + 4z - 1) 



1 - z 



ln(z) 



+ 



"ln 3 (l -z) 


+ 12 


'ln 2 (l - z)~ 




"ln(l - z)" 




1 - z 


+ 








+ 



s /241 7tt 2 , , N 
2 j i + ^- Z K^ + ^" 27rln(2) 

-|4+^-16C(3)) 



5tt 4 13C(3) \ 
36 2 / 



1 



1 + z 5 



+ — I 8Li 3 (-z) - 4Li 2 (-z) ln(z) + \ ln(l - z) ln 2 (z) + ln(z) 

1 — z \ 2 o 



2(1 + 7z 2 ) 
1 - z 



Li 3 (z) 



7 - z 2 8z 2 

• (Li 3 (l - z) + Li 2 (z) ln(l - z)) - Li 2 (z) ln(z) 



9z 2 + 7 , 3 17z 2 + 25 2 



1 - z 



1 - z 



12(1 - z) 



ln 3 (z) 



2(1 -z) 



1 - z 

4C(3) 
1 - z 

ln 2 (l-z)ln(z) 



x1 >, s n , 8z(z + l) T , 21z 2 + 22z-22 T , 33z 4 - 40z 3 + 3z 2 - 18z + 6 2 , 
4(z + 1) In (1-z) Li 2 (-z) — Li 2 (z) H — — -= In (z) 



1 - z 



3(1 -z) 



21z 2 + 46z-13 1 , , , , 8z(z+l), , , , . . , 2 „ •>- - - », 

■ ln(l - z) ln(z) ^ ^ ln(z + 1) ln(z) - 2(z + 6) ln 2 (l - z) - — -tt 2 ln(l - z) 



12(1 -z) 3 

3z 2 -5 



3(1 -z) "~ v ~ -,- v -, x _ z 
63z 4 - 144z 3 + 40z 2 - 32z + 9 , , . Voz - 5 , . 
6(l-z)»(z + l) ln(z) + —2~ ln(1 " Z) + 



2(1 -z) 

3z 2 -10z + 6 2 21z 2 -26z + 45 



6(1 -z) 



-TT- + ■ 



6(1 -z) 



(60) 



Q 



U r/ n x 11 



1 - z 



22 


'ln(l-z)" 


+ 


/34 7T 2 \ 


1 


1 + z 2 


y 


1 - z 


+ 


U 3 J 


1 - z_ 


+ + 2(l-z) 



+ + { 4(1 -,,(» + l£ -3 C ,3, 



' ln2(z) + 6(T^zy ln(z) + y {1 + z) ln(1 ~ z) 



7T 



+ + 



9 J 

68 2+^ [ln(l - z " 
ln 3 (z) 



43 - 77zl T 

J>i + (5(l-z) 



/ H41 41„ 2 
V 288 54 



., tt 4 C(3)\ 22 

+ 77 . 2 + 7 r 2 ln(2 ) + --^i + y 



ln 2 (l-z) 



1 - z 



1 - z 



55 tt 2 

27 + y- 9C(3) 



1 



1 - z 



1 + z 2 

1 - z 



(-5Li 3 (z)-4Li 3 (-z)-Li 3 (l-z) 



+ ^^ + Li 2 (z)(21n(z) + ln(l - z)) + 2Li 2 (-z) ln(z) - ln(l - z) h+(z) + ln 2 (l - z) ln(z) + — ln(z) 

5z 2 -13^ oN tt 2 (z 2 -3) 1 . . 33z 2 -2z + 8 T . . , Az(l + z) T ., . 2z 2 + 4z + 7 , , 

-C( 3 ) + ~Th r ln ( X - z ) + - * — Li2 ( z ) + - . -Li 2 (-«) ^rq r- ln(l - *) hi(z) 



2(1 -z) 



6(1 -z) 



6(1 -z) 

13z 4 -18z 3 + 2z 2 -2z-11 2 , , 4z(z + l), (412 + 47). a,, , 
ln 2 (z) H ^ ln(z + 1) ln(z) - ln 2 (l - z) 



1 - z 



24(1 - z) 3 



1 - z 



3(1 -z) 

6z 2 - 4z + 3 



12 



12(1 - z) 



173z-37, ,„ , 43z 4 -159z 3 + 96z 2 -105z + 29, . , 2z 2 - 39z - 53 
+ — ln(l-z) + — ln(z) + 



18 



18(l-z) 2 (z + l) 



27(1 - z) 



(61) 



12 



Q 




1-z 



+ (l + z)ln(z) 



00 i 




8 


"ln(l - z) 


8 


1 


I V2 


f 3 


1-z 


9 


1 - z_ 



Az z + z - 5 



+ (l + z)ln 2 (z) + 




(8z 3 + 17z 2 - 12z - 17) 



1, Nl . . 56z 2 + 52z-80 201 p . , /3139 tt 2 2 
+ 1) ln(l -z) + z — >L + 5(l-z)[ 



9,l i + 5(1 -^lw-T + l c(3) 



3" ' K ' 9 
+ (z + 1) f-8Lis(«) - 16Li 3 (-z) + 41n(z)(2Li 2 (-z) + Li 2 (z)) + 
4(z - l)(z 2 +4z + 1) 



3(1 


-z) 


\ 56 


1 


j + 27 


1 - 2 . 



+ 



+ 



ln 3 (z) 
6 



- 4C(3) 



3 z 



(2Li 2 (-z) + Li 2 (z) + ln(l - z) ln(z) + 2 ln(z + 1) ln(z)) + 



20z 9 + 45z 7 + 40z 8 - 339z 6 - 215z 5 - 345z 2 + 921z 4 - 457z 3 + 159z + 43 

9(z + 1) 3 (1 -z) 4 

400z 7 + 93z 6 - 1970z 5 + 963z 4 + 1396z 3 - 765z 2 - 11702 + 669 56 
54(z + 1) 2 (1 -z) 3 + 27? 



ln(z) 



ln(z) 



iz 3 + 17z 2 - 12z - 17 
12(1 - z) 



ln 2 (z) 



(62) 



F, 



(C F T R ni 



Q 



1 - z 



1 - z 



1 + z 



2(1 4- z 2 ) 4 

^-iln(z)--(l + z)ln(l-z) + 




"ln 2 (l-z)~ 


16 


"ln(l - z)~ 


4 


1 


1 + z 2 / 


1 - z 


+ a 


1 - z 




1 - z_ 


( 



-ln(l-z)ln(z) 

4 / ,m a,i x 5z 2 -12z + 5 w , 16 ,„ 1X1 „ , 19z-23 
+ - (z + 1) ln 2 (l - z) 9(1 _ z ) ln(z) - Y { ~ } ( ~ ] + 27 " 



-*)H 


8 
"3 


'173 


8 


,^2~~ 


"27 



ln(l - z) 



1-z 



ln 2 (z) 



(63) 



(2) (2) 

The results for the functions cLq and a^' presented above satisfy the fermion number conservation condition 
Ea. (|36|l . The two functions are not separately integrable over the interval < z < 1 because of the terms ~ 1/z; 
however, all such terms cancel in the difference. The integral of the function p^ FTRni "> ovcr < z < 1 vanishes as 

(2) 

expected, since d— has no terms proportional to ni . 

Our result can be compared with the large-/?o approximation of Ref. |2^. We have checked that the 0(/3oa 2 ) term 
in that paper coincides with the n^-dependent part of our result Eci. (|59"| . Another check of our result is provided 
by the soft, z — > 1, limit. In such kinematic regime, almost all the energy of the initial parton is transferred to the 
observed heavy quark in the final state; the energy radiated away in the fragmentation process is small. This leads 
to large contributions to D ini , of the form a™ln fc (l — z)/(l — z) with k < 2n — 1, which can be resummed using 
the soft gluon resummation formalism [||. When the resummed result in~0| is expanded in powers of a s through 
(D(a 2 ), an approximation for S- 2 \ valid for z rs 1, is obtained. We can compare those predictions with the explicit 
calculation presented in this paper. It is convenient to perform such comparison considering the Mellin moments of 
S- 2 \z). Recall that the Mellin transform is defined as: 



4f(iV) = dz z N - 1 d i Q ) {z). (64) 
Jo 

A useful summary of Mellin moments is given in Ref. 24] . The soft, z — > 1 limit corresponds to N — > oo limit in the 
Mellin space. In what follows, we present the non-vanishing asymptotics of d Q (N), for N — ► oo: 

2 , 



F, 



V») -,21n 4 (A0 + (4L + 8 7E 



4) In (TV) 



2L' 



(12 7E -7)L+-7r- 



12 7 2 - 12 7E ) In 2 (TV) 
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+ ((4 7 e - 3)L 2 + \^ir 2 - 1 + 12 7 2 - 14 7E j L + ^(tt 2 - 6 + 67I - <r„ )(2-. L - 1 ) ) ln( V i 
+ ^(47E - 3) 2 L 2 + (-tt 2 + y + 6C(3) - 7 E - 7 7 1 + ^p- + 4r^ j I + J-. L - - 2;r ln(2) - H 
- §C(3) + ^ + 2 7 | + V 7 | - V 7E - 2 7 2 - ^ + ™, 

F^ CF) - -^ln 3 W + (-%L - 2 4 7E - 3 4 + 4) ln 2 (A) + ■ ' 34 ■ ^ 



9 V/ V 3 3 ' 9 3/ v/ V 6 V 9 3 3 
55 14 2 n>/ . 68 2 2 22 ,\ , AT , / 11 7E 11\ r2 A 1:L 7f 3 5 

" 27 " <T + 9C(3) - <T E + r ^ - T 7e j ln(7V) + (~f + T J L + + T " 3C(3) 



9 3 J 9 288 9 54 3 12 9 18 27 

^> - I ^ + (?* + 5,3^ + !)!-») + - i) * + + £ - §) , 

_56 7E 3139 _^ 2 2 
27 + 648 3 + 3 U j ' 

§ + (ix + § + f ™) i»» + 0i 2 + + 1) £ + M 7E _ 1 + 1 7| + j^) ,„ (N) 

Here 7E is the Euler constant. 

Comparing these results with Ref . |6j , we find full agreement for all terms a s (Ho)ln k (N),k> 2 that were investigated 
there, provided that we use the matching relation for the coupling constants evolving with ni and ni + 1 flavors. The 
results for the subleading ln fc (A), k — 0, 1 terms are new; they can be used to extract the coefficient (cf. Eq.(70) 
in 6] ) needed to extend the soft gluon resummation for D ml to next-to-next-to- leading logarithmic accuracy. 



VI. CONCLUSIONS 



In this paper, we derived the initial condition for the heavy quark perturbative fragmentation function through 
0(a 2 ). The importance of this result is twofold. First, it enables us to derive energy distributions of heavy quarks 
produced in hard scattering processes from pure massless calculations. In addition, it can be used to resum large 
collinear logarithms through NNLL level, using the DGLAP evolution equation. 

To derive the initial condition, we make use of the process-independent approach that was recently proposed in 
0, QJ ■ We suggest a simple modification of the original proposal that renders the methods for multi-loop calculations 
applicable to the calculation of the initial condition of the perturbative fragmentation function. 

In the context of QED, an interesting potential application of the initial condition, derived in this paper, is the 
calculation of the electron energy spectrum in muon decay. Currently, this observable is being very accurately 
measured in TWIST experiment. The goal of the experiment is to extract Michel parameters |2^| from the electron 
energy spectrum with a relative precision at the level of 10 -4 . Known results on 0(a 2 ln™(m^/m e )), n — 1, 2 terms 
0, |5g suggest that 0(a 2 ) corrections without logarithmic enhancement might be important for an unambiguous 
interpretation of the experimental result. With the initial condition for the fragmentation function at hand, one 
can obtain the electron energy spectrum in muon decay by computing the energy spectrum for massless electron, 
removing collinear divergencies by MS renormalization and convoluting decay rate, obtained in this way, with the 
initial condition for the fragmentation function derived in this paper. This approach simplifies the calculation of the 
electron energy spectrum in muon decay. Similar methods can be used to extend the results of Ref. 27] on QED 
corrections to deep-inelastic scattering. 

In the context of collider physics, there are many applications for the initial condition of the perturbative heavy 
quark fragmentation function. First, the B-meson energy spectrum was measured by ALEPH, OPAL and SLD 
collaborations; this is currently the primary source of information about the 6-quark fragmentation function. It 
has been observed that inclusion of higher order QCD corrections reduces the significance of non-perturbative 
effects in the bottom quark fragmentation. The knowledge of D lm through 0(a 2 s ) : permits a reanalysis of B- meson 
fragmentation with the NNLL accuracy. Another example of potential application is an accurate prediction of the 



14 



6-quark energy spectrum in top quark decay. This is currently known to 0(a s ) |29| . The calculation presented in this 
paper offers the possibility to extend this analysis to 0(a 2 s ). 
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